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Abstract 



We present a high statistic systematic study of the overlap correlation function well 
below the critical temperature in the three dimensional Gaussian spin glass. The off- 
^ ' equilibrium correlation function has been studied confirming the power law behavior 

for the dynamical correlation length. In particular we have computed the dynamical 
critical exponent z in a wide range of temperatures, 0.35 < T < 0.9, obtaining a de- 
pendence z{T) = 6.2/T in a very good agreement with recent experiments. Moreover, 
we report a study of the violation of the fluctuation-dissipation theorem for very low 
temperatures T = 0.5 and T = 0.35. All our numerical results avoid a droplet model 
interpretation even when T is as low as T = 0.35. 



1 Introduction 

The nature of the low temperature phase of finite dimensional spin glasses is still a subject 
of controversy i a ill. 

Recently Bray, Moore, Bokil and Drossel |^, ^ have questioned many of the numerical 
results obtained with Monte Carlo methods in the three dimensional Edwards-Anderson 
(EA) model i,0,ii0. 

Inspired by the study of the Migdal-Kadanoff approximation (MKA) of the EA model, 
they argued that the numerical results, that were obtained at temperatures T > 3/4 Tc, 
could be strongly affected by finite size effects and that one should go to sizes larger than 
the crossover length L* in order to see the right (droplet) behavior. They found (in the 
framework of the MKA) that the crossover length is L* ~ 100 for T ~ 0.7 Tc and that it 
decreases for lower temperatures: L* < 10 when T < 0.5 T^ P, |^ (see also the comment 

This is maybe the main motivation that pushed us to study the EA model in the very 
low temperature region: verify whether the behavior already found at T ~ 0.75 Tc persists at 
T < 0.5 Tc. In fact at these temperatures we can simulate (using off-equilibrium techniques) 
a system of size larger than L* (in this paper, we will present data for sizes L = 24 and 
L = 64). 

The numerical data showed in this paper have been measured in the off-equilibrium 
dynamical regime. This way of probing the system properties, apart from being much more 
similar to the experimental procedure, does not present the thermalization problems of a 
simulation performed at the equilibrium, that would be insurmountable obstacles at so low 
temperatures. The efficiency of this way of measuring has been largely tested in the recent 



past 0, |1^, 0, |Tg, |I^. Moreover the off-equilibrium dynamics of spin glasses has received in 
the last years a great attention both from the experimental |p!7| , |18| and from the analytical 
19| points of view. 



Taking the measurements in the off-equilibrium regime we are able to confront the droplet 
model (DM) [|l| and the Mean Field like theory (MP) 0] on two grounds: the off-equilibrium 
regime itself and the equilibrium one, that can be obtained in the limit of very large times. 
We can take this limit quite safely thanks to the very large time reached in our simulations. 

A preliminary analysis based on the data at temperatures T = 0.7 and T = 0.35 was 
reported in reference |T^. In this paper we present an extended analysis based on nine 



different temperatures obtaining a precise temperature dependence of the dynamical critical 
exponent in order to have an accurate confront with recent experiments. 

2 The model and the numerical method 

We have simulated the Gaussian Ising spin glass on a three-dimensional cubic lattice of 
volume L^ with periodic boundary conditions. The Hamiltonian of the system is 



'H = - Yl ^»^' 



ij'^j 



<«j> 



We denote by < ij > the sum over nearest neighbor pairs. Jij are Gaussian variables with 
zero mean and unit variance. 



We focus our attention on the study of the point-point overlap correlation function com- 
puted at distance x and time t 

G{x, t) = —Yl {cri+xri+xCriTi)t . (2) 

where a and r are two real replicas (systems which evolve with the same disorder) and the 
index i runs over all the points of the lattice. As usual we denote by (■ ■ ■) the average over 
the disorder and, in this context, ((■ ■ ■))j is the average over the dynamical process until time 
t (for a given realization of the disorder). The two replicas {a and r) evolve with different 
random numbers. 

The simulation has been performed in a similar way to the experimental procedure: the 
system is prepared in a high temperature configuration (actually the initial configurations 
were chosen at random, i.e. T = cxd) and suddenly it is quenched below the (estimated) 
critical temperature, T^ = 0.95(3) ^. Immediately we start taking the measurements, 
which obviously depend on time. The equilibrium behavior is recovered in the large time 
limit. We have used as dynamical process the standard Metropolis method. 

We have simulated 4 samples (8 systems) of an L = 64 lattice, measuring the correlation 
function at times t = 100 ■ 2^^ (with k = 0, . . . , 13) and temperatures T = 0.9, 0.8, 0.7, 
0.6, 0.5, 0.4, and 0.35. In addition, we have simulated 4096 samples of an L = 24 lattice 
measuring at times t = 2^ (with k = 7, . . . , 19) and at three temperatures: T = 0.7, 0.5 and 
0.35. 

For the study of the fluctuation-dissipation relation we have used L = 64 systems and we 
have simulated them for more than 10^ Monte Carlo steps. All the simulations have been 
performed with the help of the parallel computer APEIOO [pO[] . 

3 Results on the correlation function 

First, we analyze the correlation functions computed with the L = 64 runs. The study of 
the numerical data suggest us the following Ansatz for the time and spatial dependences of 



the correlation function 12 



const 
G{x,t) = — — exp 



X \ 



Xit), 



(3) 



where ^(t) is the dynamical correlation length. The numerical data clearly show that the 
dynamical correlation length depends on the time following a power law ^(t) = Bt^^^ where 
z is the dynamical critical exponent. The exponents a, 5 and z and the amplitude B could, 
in principle, depend on the temperature. However we obtain (see below) that a, 5 and B 
are almost temperature independent, while z{T) is inversely proportional to T. 

In Table ^we report the results of our fits (always done using the CERN routine MINUIT 
||21||). We remark that, for a given temperature, we have fitted our numerical data to the 
Ansatz of Eq.(^) in two steps. In the first step we fix the distance in the correlation function 
and we perform the following three parameters fit in the variable t 

log G{x, t) = A{x) - B{x) r^/" . (4) 



T 


z{T) 


6 


B-' 


a 


0.9 


6.85(1.0) 


1.37(11) 


1.02(4) 


0.60(7) 


0.8 


7.5(1.3) 


1.37(9) 


1.06(4) 


0.49(9) 


0.7 


9.3(0.7) 


1.50(5) 


1.02(4) 


0.53(5) 


0.6 


10.3(1.2) 


1.38(3) 


1.09(4) 


0.49(16) 


0.5 


11.7(1.8) 


1.43(2) 


1.04(4) 


0.59(20) 


0.4 


14.1(2.4) 


1.45(4) 


0.99(3) 


0.60(26) 


0.35 


19.9(3.8) 


1.41(6) 


1.03(7) 


0.29(32) 



Table 1: Parameters of the point-point correlation function 



We have found that 6/z is independent of x. The second step has been to extract from 
A{x) and B{x) the exponents a and 6 and the amplitude B using the formulae: A{x) = 
const — alogx and B{x) = B~^x^. We report our final values of z, 6, B and a in Table 0. 

The resulting values for z{T) (see Table 0) can be fitted to a power law (using all the 
temperatures of Table ^ obtaining 



z(T)= 6.4(6) T-°-9^(20) 



(5) 



From the previous fit we can guess a simpler law for the dynamical critical exponent z{T) = 
a/T, obtaining^ 

<T) ^ ^ . (6) 

This kind of behavior suggests that the low temperature dynamics in spin glasses is domi- 
nated mainly by activated processes with free energy barriers diverging logarithmically with 
the size of the system. 

We can finally write down the dependence of the dynamical correlation length on the 
time as well as on the temperature: 

e(t, T) OC t^/^-2(3) = t0161(8) T ^ ^0.153(12) T/T. ^ (7) 

where we have assumed that the temperature of the phase transition is Tc = 0.95(3) |^. 
The agreement of the previous formula with the experiments is very good. We recall that in 
experiments |T7] it was found the following dependence for the dynamical correlation length 

e(t,r)oct°-i^°^/^'' . (8) 

where Tg is the experimental critical temperature (the authors of this result do not quote 
the error in the exponent). 

A further check of Eq. (^ would be the collapse of the data (measured at different times 
and different temperatures) when plotting G{x,t)x°' versus x/t^^^^'^\ To this purpose we 
use the data from 4 samples of the 64'^ runs, together with those measured on 4096 samples 
of 24^ runs. We remark that, in the 24^ runs, the volume is near 19 times less than the 
L = 64 runs but we have computed 1000 times more samples and so we expect the errors to 
be smaller. 



^ This law was found by Kisker et al. for the ±1 three dimensfonal spin glass (See references [^2|, |2^). 
Moreover this law was guessed for the Gaussian model using numerical data taken at temperatures T — 0.7 
and r = 0.35 in reference ||l31. 
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Figure 1: Scaling plot for the correlation function G{x,t) measured at two very low temper- 
atures, T = 0.5 and 0.35, and two lattice sizes, L = 24 and 64. It shows that the finite size 
effects are negligible and it also gives reliability to our estimate for ^(t). 



L=24 







Figure 2: The scaled data for G(x, t) are well described by a temperature-independent scaling 
function, that can be very well approximated by our fit. 



In Fig. |1] we plot the correlation function for two low temperatures (T = 0.5 and 0.35) 
using as variables x/^{t) and x°^G{x,t) (we have taken a = 0.5, see Table jT]). In the plots 
we use the data from both runs {L = 24 and L = 64) and they superimpose perfectly. In 
the insets we present the same data in a log-linear scale in order to let the reader evaluate 
better the collapse. It is clear that the scaling is impressive even at the lowest temperature 
T = 0.35. We can also state that the finite size effects are negligible for the lattice sizes 
used. 

Scaling arguments tell us that the more general scaling function for the correlation func- 
tion is (for large x and t) 

G{x,t)^x-"g(^-^^ , (9) 

where the scaling function g{y) is smooth. Moreover, in the scaling regime, g{y) should not 
depend neither on temperature nor on the lattice size. Note that in our Ansatz (|^) we have 
chosen an exponential function for the scaling function: g{y) oc exp{—y^), and we show that 
it fits very well the data. However to check that our estimates of a and C,{t) are correct we 
do not need to know g{y)- We can simply plot x"G(x,t) versus x/^(t) (as it was done in 
Fig. |I]) and check how well the data collapse. 

In order to check the temperature independence of g{y) we show in Fig. ^ the scaling 
function for three different temperatures (T = 0.35,0.5 and 0.7), together with the expo- 



nential function exp[— y 



1.42(2)1 



(see Table |l]) obtained through the fitting procedure. It is 



clear that the scahng function is really temperature-independent and it can be very well 
approximated by the exponential function as we have chosen in our Ansatz. 

Another interesting issue is the extrapolation of the correlation function to infinite time. 
In this limit we can compare again our numerical results with the predictions of the droplet 
model and with that of the RSB theory. In the former the extrapolated correlation function 
tends to the value g^^ for large distances, whereas the RSB prediction is a pure power law 
going asymptotically to zero p4|. Our Ansatz, which describes perfectly the numerical data. 



supports the RSB prediction even for the lowest temperatures. 

Nevertheless, we have tried to fit our data with a functional dependence compatible with 
the droplet model, that is, G{x,t) = Gooix) g{x/^{t)), where Gooix) = Ax'^' + G. li G = 0, 
then the previous formula is exactly our Ansatz (and it implies a breaking of the replica 
symmetry), while if C = g^^ then it would support a droplet picture. Fitting the data to the 
previous formula, Goo{x), we have found that at every temperature and even at T = 0.35 
(our lowest temperature), the best value for G is always compatible with zero. At very 
low temperatures, i.e. T = 0.35, the Edwards-Anderson order parameter is so close to one 
(g^A — 1) than we can safely distinguish between the two competing theories. In fact in 
the droplet like formula we have that God{x) is almost constant^ and so we simply should 
fit the data into the scaling formula of Eq.(|^) without the factor a;~" in order to check the 
correctness of the droplet model. 

In Fig. 1^ we present the data rescaled with the formula suggested by the droplet model 
(left plot) and with that implied by RSB (right plot). It is clear that the RSB prediction 
fits much better the numerical data. Note that the data error is sufficiently small to affirm 
safely that the data in the left plot have no collapse at all. A scaling plot like the left one 



has been recently presented by Komori et al. in [^ (see also reference [^). We believe that 
the rather poor collapse of their data (see Fig. 5 in [|^) is due to the fact that they neglect 
the factor a;~" in the scaling formula. A much better collapse would be obtained by plotting 

y/x G{x,t) versus x/^(t) (see p6[| ). 



4 Fluctuation dissipation relation at very low temper- 
atures 

Now, we present the results of the analysis based on the generalization of the fluctuation- 
dissipation theorem (FDT) in the out of equilibrium regime ||2^ . In this section we will focus 



on the scaling properties of the aging region and the violation of fluctuation-dissipation at 
very low temperatures. 



A preliminary analysis was done in [28|| studying the violation of FDT at temperature, 
T = 0.7. Here we have simulated different lower temperatures and so, as byproduct, we 
can study the scaling properties of the violation of FDT. An analogous analysis, with many 



temperatures, was done in [28| but on the four- dimensional EA model. 



For the sake of conciseness, we do not repeat all the formalism and we address the 



interested reader to one of the previous publications on the subject [^ |2^, ^ |30|, |T^. Here 



^ Actually it slowly decreases from 1 to g^^ ~ 1, but for all our purposes it can be considered as a 
constant. 



10' 



10' 



10' 



10' 




Figure 3: The comparison of these two scahng plots clearly shows that the -/x factor is 
essential in order to collapse the correlation function data. The temperature is T = 0.35 and 
the lattice size L = 24. 



we simply recall the main formulae that we use. As usual we define the integrated response 
to a very small external field as 

X{t, t^) = hm -^ /* R{t, t') h{t') dt' , (10) 

where h{t) = h^ 9{t — t^) and R{t, t') = -^J2i Qhu/) ■ The autocorrelation function is defined 

as 

C{t,t^) = ^J2i^^it)s,{t^)) . (11) 

i 

Relating these two functions, in the large times limit, via 

Tx{t,t^) = S{C{t,t^)) , (12) 

we have that, at the equilibrium, the fiuctuation-dissipation theorem (FDT) holds and 
S{C) = 1 — C, while in the aging regime the function S{C) can be linked to the equilibrium 



overlap distribution through P{q) = g^ 



Models that, in the frozen phase, do not show any breaking of the replica symmetry, 
have, at the equilibrium level, a static P{q) = S{q — ^ea), which dynamically corresponds 
to the absence of response in the aging regime. This means that, plotting xit^tw) versus 
Cit^tw), we obtain a horizontal line in the range C < q^A |23 (in the quasi-equilibrium 



regime, C > q^A, and it always holds Tx = 1 — C independently of the model) 

In Fig. 1^ we show the results for different temperatures in the usual plot x(t, t^) versus 
C{t,tuj). Note that in this plot the FDT line is x = (1 ~ C')/T and so it is different for 
different temperatures. It is quite clear that, even for very large times, the curves are far 
from been horizontal when they leave the FDT line. This result gives more evidence in favor 



of a replica symmetry breaking in the very low temperature phase of the 3D EA model [^ . 

We present the data for different temperatures on a single plot in order to make more 

evident the fact that the numerical data seem to stay on the same curve once the system 

enters into the aging regime, i.e. when the points leave the FDT line. This kind of behavior 



has been observed in the four-dimensional EA model |^ and it is reminiscent of the mean- 
field solution. 

Indeed in the SK model, using the Par isi- Toulouse (PAT) hypothesis P^, it can be shown 
H that 

Q(n\ - I ^~'^ forC>gEA(T) , . - 

^^^>-\TVT^~C ioiC<q,A{T) . ^ ' 

The formula can be easily generalized assuming a generic power law behavior in the aging 
regime: S{C) = T A{1 — C)^ (the mean-field value for the exponent is B = 1/2). 

We use this generalization to fit the data and we obtain very good results. The best fit 
parameters have been estimated from the collapse of the data reported in Fig. ^ and they 
are A ~ 0.7 and B ^ 0.41 (to be compared with the mean-field values A = 1 and B = 1/2, 
and those obtained for the 4D EA model A ~ 0.52 and B ~ 0.41 ^). 

In order to show the validity of the fitting formula, we present in Fig. ^ the collapse of 
the scaled data using the variables a; = (1 — C)T~^ and y = xT^~'^y where (p = j^ = 1.7. 
It is easy to see that, if the previous scaling holds, the data should stay on two power laws: 
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Figure 4: Response against the autocorrelation function for three different temperatures and 
lattice size L = 64. Note that in this figure we plot x(^, tw) versus C(t, tw). The data stay on 
a single universal curve when they leave the FDT lines. This curve is clearly non horizontal 
and this hints for a breaking of the replica symmetry in the very low temperature phase of 
the EA model. 
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Figure 5: Scaling plot of the curves shown in the previous figure: xT'^~'^ versus (1 — C) T'"^ 
with = 1.7. In the FDT regime (left part of the figure) the scaled data stay on the line 
y = X, while in the off-equilibrium regime (right part of the figure) they follow the plotted 
power law y = Ax^ with B = 0.41. 
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y = X and y = Ax^ in the quasi-equilibrium and aging regime, respectively. The two power 
laws are reported in Fig. ^. 

Even if we may expect a breakdown of the assumed scaling for large values of the scaling 
variable x (i.e. the scaled data are no longer described by a power law), we note however that 
for a quite large range the collapse is very good and very well approximated by a power law. 
Moreover, we remark that the collapse has been obtained adjusting only one parameter. 

5 Discussion 

We have studied the off equilibrium dynamics of the three dimensional Gaussian spin glass 
in the very low temperature phase. In particular we have studied the scaling properties of 
the dynamical overlap correlation functions and the scaling properties of the violation of the 
fluctuation-dissipation. 

We have tried to fit our correlation functions to the functional form predicted by the 
droplet model but the fits were poor. Moreover a correlation length diverging following a 
power law with the time implies, as was noted by Rieger [^], barriers diverging not as L^ 



(as predicted by the droplet model with the lower bound ijj > 6 c::^ 0.2) but as logL. This 
latter results implies ip = 0, hence violating the droplet lower bound. 

It is interesting to note that the experimental data could be fitted to the droplet formula 



assuming that ip = 9 |T^. However while both, the results of numerical simulations and the 
experiments are in very good agreement with a power law fit for ^{t,T), the numerical fit 
assuming a droplet formula for ^(t,T) |Q disagrees with the experimental fit assuming the 
same hypothesis |]T^ . 

As it has been noted above, our final result for the dynamical correlation length is in a 
very good agreement with the experimental result. 

We remark that the same scenario (power law dependence of ^(t, T) and linear dependence 
of 1/z with temperature) also emerges in four and six dimensions. In the latter case it was 
found that z{T) = AT^/T |I^ (z = 4 at the transition is the value predicted by Mean 
Field) while in the former one z{T) = 5.5 Tc/T |]13[. Moreover in these two dimensions the 



overlap correlation function constrained to zero overlap follows a pure power law as in three 
dimensions. 

If we send to infinity the time in our Ansatz for the overlap-overlap correlation function 
we obtain a pure power decay G{x) oc x~" with a ~ 0.5, with a small dependence of a on the 
temperature for the whole spin glass phase. We recall again that the droplet prediction is 
G{x) — >■ q^^ in contradiction with our numerical correlation functions (this fact was already 
noted in [0). Instead, the pure power behavior is supported by the Gaussian approximation 



using the Mean Field solution |24 



One could argue that the simulated temperatures are not low enough and the times 
and sizes not large enough in order to see the "true" (droplet) behavior of the EA model. 
However, as we stressed in the Introduction, L = 64 is large enough for temperatures as low 
as T = 0.35 and T = 0.5. Moreover our large times extrapolations are very safe thanks to 
the measurements have been taken over 6 time decades. 

We have shown numerical results that contradict the droplet predictions in a wide range 
of temperatures (0.35 < T < 0.9). In particular we point out that our results (both for 



12 



correlation functions and for violation of FDT) at a very low temperature, T = 0.35, support 
a Mean Field picture. 

Finally, we remark that using the (PAT) Mean Field scaling relations for the P{q) |^ we 
have obtained a very good scaling plot of the violation of fluctuation-dissipation (like in four 



dimensions p8|). This provides us another strong evidence calling for a low temperature 



phase being well described by Mean Field [Q. 
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